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Abstract 

We consti uct static and mass-shedding limit sequences of hybrid stars, composed of colour flavour locked 
(CFL) quark matter core, for a set of equations of state (EOSs). The EOS for the hadronic matter is obtained 
using appropriately calibrated extended field theoretical based relativistic mean-field model. The MIT bag 
model is employed to compute the EOSs of the CFL quark matter for difi"erent values of the CFL gap 
parameter in the range of 50-150 MeV with the deconfinement phase transition density ranging from 4po - 
6po (po - 0.16 fm~^). We find, depending on the values of the CFL gap parameter and the deconfinement 
phase transition density, the sequences of stable configurations of hybrid stars either form third families 
of the compact stars or bifurcate from the hadronic sequence. The hybrid stars have masses 1.0 - 2.\Mq 
with radii 9 - 13.5 km. The maximum values of mass shedding limit frequency for such hybrid stars are 
1-2 kHz. For the smaller values of the CFL gap parameter and the deconfinement phase transition density, 
mass-radius relationships are in harmony with those deduced by applying improved hydrogen atmosphere 
model to fit the high quality spectra from compact star X7 in the globular cluster 47 Tucanae. We observed 
for some cases that the third family of compact stars exist in the static sequence, but, disappear from the 
mass-shedding limit sequence. Our investigation suggests that the third family of compact stars in the 
mass-shedding limit sequence is more likely to appear, provided they have maximum mass in the static 
limit higher than their second family counterpart composed of pure hadronic matter. 
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I. INTRODUCTION 



Soon after the suggestion that three-flavor quark matter may be the ground state of strongly 
interacting systems 111], quark stars are postulated as possible astrophysical objects. It was also 
hypothesized that some compact stars might be hybrid stars with the core composed of quark 
matter and surrounded by nuclear mantle. The present knowledge of quantum chromodynamics 
(QCD) at high density indicates that quark matter might be in a color superconducting phases. 
The essence of color superconductivity is quark-quark color superconductor and driven 

by Bardeen, Cooper and Schrieff"er (BCS) 



3, [sl] pairing mechanism. The possible quarks color 
superconducting phases include the two-flavor color superconductor (2SC) Q, [s ] , th e colour 
flavour locked (CFL) phase and the crystalline color superconductor (CCS) 



m. 



The speculation that the colour superconducting quark matter present in the core of the hybrid 
stars has triggered many theoretical investigations. 

The hybrid stars with CFL quark matter core have been extensively studied. The hadron phase 
of the hybrid star matter is describe d by the various models which can be broa dly group ed into 
(i) non-relativistic potential models lll4ll . (ii) non-relativistic mean-field models lll5 . 



(iii) field theoretical based relativistic mean-field models (FTRMF) 



Brueckner-Hartree-Fock model 122. 



23i 
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m. 



2lh and (iv) Dirac- 



251] . The CFL quark matter appearing at the core of 
hybrid stars are described within the MIT bag model and Nambu-Jona-Lasinio (NIL) model. The 
studies based on the MIT bag model indicate the existence of stable configurations of hybrid stars 



with the CFL quark matter core 



280. Two diff'erent situations are encountered, the hybrid 



stars with CFL quark matter core either form a third family of compact stars separated from the 
purely hadronic sequence by an instability region or bifurcate from the hadronic sequence of stars 
when the central density exceeds the phase transition density at which deconfinement of hadrons to 
the CFL quark matter occurs. The scenario is completely different when NIL model is employed 
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3l|] . it was shown 



to study the hybrid stars with CFL quark matter core. Until recently [|29[ 
that the NIL like model rules out the CFL quark matter phase at the core because it renders the 
hybrid star unstable. Only very recently, it has been found that large enough values of the diquark 
coupling strength in NJL model can yield stable configurations of the hybrid star containing CFL 
quark matter core 1321 1331]. 



The stability of hybrid star with CFL quark matter core depends strongly on the values of the 
deconfinement phase transition density and the CFL gap parameter which are poorly known. In 
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the present work we construct the static sequences of hybrid stars, with CFL quark matter core, 
for a set of EOSs obtained for different values of CFL gap parameter and the deconfinement phase 
transition density. The hadron phase of the hybrid star is described by using an appropriately cali- 
brated extended FTRMF model which includes the contributions from self- and mixed interaction 
terms for cr, co and p mesons up to the quartic order. The CFL quark matter phase is described 
within the MIT bag model with an additional parameter that mimics the effect of including pertur- 
bative QCD corrections. Instead of keeping the value of the bag constant fixed as previously done 
1 2^ 27], calculations are performed for different values of the CFL gap parameter at fixed values 



of the deconfinement phase transition density. This strategy should enable us to asses better the 
influence of the CFL gap parameter on the properties of hybrid stars with CFL quark matter core. 
The CFL gap parameter A is varied in the range of 50 - 150 MeV by keeping the deconfinement 
phase transition density p, fixed in between 4po - 6po (po = 0.1 6fm"^). For the different values 
of the CFL gap parameter considered, the average quark chemical potential at the deconfinement 
phase transition density lie in the range of 375 - 500 MeV which is in reasonable agreement with 
the predictions of the NIL model. 

The paper is organized as follows, in Sec. II we describe in brief the models used to construct 
the EOSs for hadronic phase, CFL quark phase, and the mixed phase. In Sec. Ill we present 
the results for static and mass-shedding limit sequences for hybrid stars. In Sec. V we state our 
conclusions. 



II. EQUATIONS OF STATE FOR HYBRID STAR MATTER 

We construct the EOS for the hybrid star matter which is composed of hadrons at low densities, 
quark matter in the CFL phase at high densities and the mixed phase at the intermediate densities. 
The EOS for the hadron matter is obtained within the framework of the extended FTRMF model. 
The EOS for the quark matter in the CFL phase is obtained using the MIT bag model. The EOS 
for the mixed phase is constructed using the Gibbs conditions. For the hadron matter at very low 



densities p ~ 0.5po fm going down to p = 6.0 x 10 fm we use Negele-Vautherin [|34ll and 
Baym-Pethick-Sutherland EOS 
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A. Hadron phase 



The hadronic phase is described using the extended FTRMF model which includes the con- 
tributions from self- and mixed interaction terms for cr, oj, and p mesons up to the quartic order. 
The mixed interaction terms involving the p-meson field enables one to vary the density depen- 
dence of the symmetry energy coefficient and neutron skin thickness in heavy nuclei over the wide 



range without affecting the other properties of the finite nuclei [j36 



3711 ■ The contribution from 



the self-interaction of a»-mesons plays an important role in determining the high density behavior 



of EOS and consequently the structure properties of compact stars [|38L 3^. The contributions of 
self-interaction of p-meson are ignored as they affect the ground state properties of heavy nuclei 
and compact stars only very marginally [|39ll . In our recent work 113 811 we have obtained several 
parameterizations of the extended FTRMF model in such a way that the bulk nuclear observables 
and nuclear matter incompressibility coefficient are fitted well. These different parameterizations 
produce different behavior for the EOS at high densities. 

The energy density of the hadron phase in the extended FTRMF models is given by 

6Hp(Mn,l^e) = — VF + M*^dk + gojNCOiPp +pn) + ^gpNpipp - Pn) + IZmlcT^ 
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The pressure of the hadron phase matter is given by 
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where, M* = M - go-NC is the effective mass of nucleon with M being the free nucleon mass. In 
Eqs. ^ and Q, cr, co and p represent the meson fields. The g^-N, gaiN and g^^ are the meson- 
nucleon coupling strengths. The nia-, and lUp are the masses for the cr, to and p mesons. The 



coupling strengths for the self-interaction terms for the cr and to mesons are denoted by k, A and 
^. The constants a and a' represent the coupling strengths for various mixed interaction terms. 
The last term in Eqs. ([T]) and Q gives the contributions to the energy density and pressure from 
leptons, respectively. In Eqs. [Hand El Mn and jUe represent the chemical potentials for the neutrons 
and electrons, respectively. The chemical potentials for the protons fip and the muons //^ can be 
expressed in terms of and fi^ using the jS-equilibrium conditions, i.e., 

fin = lUp+lUe (3) 
f^e = f^i,- (4) 

Once th chemical potentials for the nucleons are known, their fermi-momenta can be obtained by 
solving the field equations for the mesons as given in Ref. [38]. The fermi-momenta k'^. for the 
leptons are obtained as, 

k'f = ^n]-m]. (5) 

In addition to the conditions of the yS-equilibrium, matter in pure hadronic phase is considered to 
be charge neutral. 



B. Quark matter in the CFL phase 



The free energy density for quark matter in the CFL phase is taken to be i40|] . 

where, ju is the average chemical potential for quarks and lUg is the electron chemical potential. 
The contribution to Eq. © from the quarks is given by. 



^quarks = ^ p^{p - fl)dp + ^ ^ + ' P^^P + ^'^^^ ~ "T" 



where, u and d quarks are assumed to be massless and s quark has the mass m^. The term pro 



portional to j/ in Eq. ([7]) corresponds to the QCD inspired correction [|4ll] . The second last term 
involving the CFL gap parameter A is the lowest order contribution from the formation of the CFL 
condensate. The last term B is the bag model constant which accounts for the energy difference 
between the perturbative vacuum and the true vacuum. The number densities for all the three 
flavours of quarks considered are the same and can be obtained as, 

p, = l(y^ - c^i' + 2^^lx) (8) 



with q = u,d and s and the common fermi momentum v given as, 



O / 2 ^ 



(9) 



The contribution to the Eq. ^ from the Goldstone bosons arising due to the breaking of chiral 



symmetry in the CFL phase is evaluated using fl^pj^ip., jUg) [|42ll . 



Q^^^Ou,/z,) = --/,^Ai^ 1-^ 
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where the parameters are 



. (21 - 81n2)//2 

= ,m^- 
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(10) 

(11) 

(12) 



The total energy density and pressure for the quark phase (QP) can be calculated as, 



6qp = O-CFliii, l^e) + ^I^Pq + HeiPe + P^), (13) 

and 

Pqp = -acFL(M,Me), (14) 

It is clear from the Eq. ([8]) that the densities for the u, d and s quarks are equal to each other for 
the quark matter in the CFL phase. Thus, quark matter in the CFL phase is enforced to be charge 
neutral. The electrons are present only in the mixed phase of hadronic and the quark matter. 



C. Mixed phase 

The EOS for the mixed phase composed of hadronic and CFL quark matters is obtained using 
the Gibbs conditions. The Gibbs conditions can be expressed in terms of two independent chemical 
potentials in our case as, 

PHp(Mn,f^e) = PQp{M,IUe) (15) 

where, yu„ and fie are the two independent chemical potential with //„ being the neutron chemical 
potential. In the mixed phase the average quark chemical potential fi = )U„/3. In the mixed phase 
local charge neutrality condition is replaced by the global charge neutrality 

Xptp + (^-X)p1lp = 0. (16) 
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where x is the volume fraction occupied by quark matter in the mixed phase and p"^^ is the charge 
density. It is clear from Eq. (fT6l) that both hadron and quark matter are allowed to be charged sep- 
arately. The energy density Smp and the hadron density Pmp of the mixed phase can be calculated 



asm. 



&MP = X&QP + (1 -X)&HP, (17) 

Pmp =XPqp + {'^- X)Php, ( 1 8) 

Once these quantities are determined, we can construct the complete EOS with the hadron phase, 
quark matter phase and mixed phase and compute the properties of hybrid compact star. 

ni. STRUCTURE OF HYBRID STARS WITH CFL CORE 

We study the properties of the hybrid stars, composed of CFL quark matter core, for a set 



of EOSs obtained for different values of the CFL gap parameter A and the deconfineme nt p ' 



hase 



transition density pt. Instead of fixed values of the bag constant as customarily done [26, 
we adjust the bag constant for each values of the CFL gap parameter to yield the desired value 
of the phase transition density. We consider the values of the CFL gap parameter in the range 



33L 



45 



m 



of 50 - 150 MeV as estimated and employed for the studies of hybrid stars [|26l . 
In Fig. [Uwe plot the several EOSs for the hybrid star matter obtained for different values of the 
CFL gap parameter A with p, = Ap^ - 6po. The EOS for the hadron phase is obtained within 
the framework of the extended FTRMF model as discussed in the preceeding section. In Ref. 



13811 we have obtained several parameter sets for the extended FTRMF model for different values 
of the coupling strength of the cu-meson self-interaction term and the neutron-skin thickness in 
208p|^ nucleus as these are not well determined from the presently available experimental data. 
Each of the parameterizations are consistent with bulk properties of the finite nuclei and nuclear 
matter. In the present work we have employed the parameter set which corresponds to cu-meson 
self-interaction strength ^ = (Eq. [T]) and neutron-skin thickness of 0.2 fm in the ^"*^Pb nucleus. 
The choice of ^ = yields stiff EOS for the hadronic matter at high density. The EOSs of CFL 
quark matter for different values of A and are obtained using strange quark mass = 150 MeV 
and the constant c = 0.3 in Eq. (|7]). In Fig. [2] we plot the values of the bag constant B^^^ as 
a function of A for pt = 4po - 6po. We get the values of B somewhat higher compared to the 
ones commonly used. In particular, the value of B increases rapidly with the deconfinement phase 
transition density. This is because, stiffness of the EOS for the hadronic matter as considered 
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here increases rapidly with density. Further, as the CFL gap parameter increases, the value of B 
increases to keep the phase transition density unaltered. In Fig. [3]we plot the values of the average 
quark chemical potential fUt at the deconfinement phase transition densities pt = 4po - 6po as a 
function of the CFL gap parameter. For our choice of the phase transition densities, the values of 
yu, are in the range of 375 - 500 MeV which are in reasonable agreement with the ones obtained 



inRefs. (13, 



3211 . The properties of the static and rotating compact stars resulting from our set of 



EOSs are computed using the code developed by Stergioulas [|47ll 



A. Static sequences 

The sequence of static compact stars is obtained by varying the central energy density ec for a 
given EOS. For stable configuration, 

> (19) 

where, M is the gravitational mass of the static compact star. In Figs. |4]-[6]we plot the mass- 
radius relationships for static sequences obtained for various EOSs corresponding to the different 
values of the CFL gap parameter A = 50 - 150 MeV with p, = 4po - 6po. The solid circle on 
each of the curves marks the point at which the deconfinement phase transition from hadron to 
CFL quark matter occur. The curves on the left to the solid circles represent the sequences of the 
hybrid stars with CFL quark matter core. The black dotted line represents the static sequence of 
the compact stars composed of pure hadronic matter. The radius of hybrid stars decreases with 
increasing central energy density. It can be seen that the stable configurations of hybrid stars with 
CFL quark matter core either bifurcate from the hadronic sequence or form different branch so- 



called third family of compact stars [1431 148l1 . With p, = 4po - 5po, the stable configurations of 
the hybrid stars exist for all the values of the CFL gap parameter considered. In particular, for 
A = 50 - 100 MeV with p; = 4po, sequences of stable configurations of the hybrid stars bifurcate 
from the hadronic sequence at the central density exceeding the one at which the onset of mixed 
phase occur. For all the other cases with Pf = 4po - 5po, the hybrid stars belong to the third families 
of the compact stars. When, the value of p, is increased to 6po, the hybrid stars with CFL quark 
matter core become stable only for A ^ 125 MeV. We get the masses for such hybrid stars in the 
range of 1.0 - 2.1Mo with radii 9.3 - 13.5km. 

For the comparison, in Figs. |4]- 16] we plot contours (dot-dashed/maroon) in the mass-radius 
plane which are deduced with 90% confidence, by fitting the high quality spectra from the compact 
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star X7 in the globular cluster 47 Tucanae [|49h . These spectra were fitted within an improved 



hydrogen atmosphere model which also accounts for the variations in the surface gravity with 
mass and radius of the compact stars. These M - R contours indicate that a compact star with the 
canonical mass (1.4Mo) should have radius in the range of ~ 13 - 17 km, whereas, compact star 
with the canonical radius (10 km) has mass ~ 2Mq. The compact stars composed of only hadronic 
matter can satisfy either the constraint on the radius at the canonical mass or the constraint on the 
mass at the canonical radius li49tl . Our results for the A = 50 - 75 MeV with = 4po - 5po are 
bounded by the dot-dashed maroon contours over broad range of mass and radius. For these cases, 
compact stars with canonical mass 1.4Mo have radii 13.5km and are composed of only hadronic 
matter. But, compact stars with radii around the canonical value 10km have masses nearly 2Mq 
and are the hybrid stars with CFL quark matter core. We also plot the M - R curves obtained 
using the constraints imposed by the discoveries of the kHz quasi-periodic oscillations (QPOs) 



pOll and the X-ray transient XTE J1739-285 [51]. The frequency of the innermost stable circular 



orbit (ISCO) inferred from the QPOs limits the mass of the non-rotating compact stars to be, 

2200 Hz 

M < Mq. (20) 

^ISCO 

The values of visco are in the rang e of 1220 - 1310 Hz. The compact star radius must be smaller 



than the ISCO, which implies [150|] 



^ 19500 Hz, 

R < km. (21) 

^ISCO 

Radii limits for masses less than the upper limit scales with M^^^. The discovery of XTE J1739- 
285 suggests that it contains a compact star rotating at 1122Hz. This, imposes the constraint on 
the maximum radius of a non-rotating compact star with mass M [52], 



M 



1/3 



We see that the hybrid 



^max<9.63^ — j . (22) 

stars with M > 1.4Mo satisfy the constraint as expressed by Eq. (I22l) . It is 



5411 that mass of the compact star rotating with 1 122 Hz is equal or larger 



also found in Refs. 115 31 
than 1.4Mo. 

In Table in we give the values of maximum masses with corresponding central energy densities 
and radii for the hadronic and the hybrid stars obtained using different values of the A and p,. 
It may be noted that, for a given A and p,, the central energy density for the maximum mass 
of the hadronic star corresponds to the one at which the onset of mixed phase occurs (see also 



Fig. [T]). The values of maximum mass 1.7 - 2.1Mq for the hybrid stars are consistent with the 
currently measured maximum mass 1.76+0.20Mo of PSR J0437-4715 n55\\ obtained by the precise 
determination of the orbital inclination angle. We note, for Pf = 5po, maximum mass of the hybrid 
stars with CFL core are nearly equal to their second family counterpart composed of hadrons. The 
hybrid stars with CFL core are smaller by about 30% compared to those second family counterpart. 
Thus, the hybrid stars with CFL core are expected to rotate significantly faster in comparison to 
the hadronic stars. We also remark that our results for the mass-radius relationship are somewhat 
similar to the ones obtained using the NJL like model 132^ |331], for the hybrid stars composed of 
CFL or CCS quark matter core. 



B. Mass-shedding limit sequences 

In Fig. Ul we plot relationship between mass and the circumferential equatorial radius R^q for 
the mass-shedding limit sequences obtained for the EOSs corresponding to diff"erent values of A 
and Pf The portion of the curves left to the solid circles represent the hybrid stars with CFL core. 
The symbol cross (x) on the different curves marks the maximum mass of the hybrid star. For 
Pt = 5po with A = 50 and 75 MeV, third family of compact stars with CFL core disappears from 
the mass- shedding limit sequence, though, they exist in the static sequence as can be seen from 
Fig. [51 Similar situation is encountered for A > 125 MeV with p, = 6po (see also Fig. (6]). In Table 
nil we give the values of the central energy density, radius and Kepler (mass-shedding) frequency 
/k at the maximum mass for the hadronic and hybrid stars corresponding to different A and p,. 
The values of the Keplerian frequency at the maximum hybrid star mass is in the range of 1 .7 - 2 
kHz, whereas, the Keplerian frequency at the maximum mass for the hadronic stars is ~ 1 kHz. 

We now examine several cases for which third family of compact stars with CFL core appears 
in the static sequence, but, disappear from the mass-shedding limit sequence. We observe that for 
these cases, maximum mass of the third family compact stars is lower than their second family 
counterpart. In Fig. [H we have plotted the mass-shedding limit sequences (upper panel) and the 
static sequences (lower panel) for A = 150 MeV with p, ranging from 5po - 6po. We clearly see 
that the third family of compact stars tend to disappear beyond Pf > 5.5po. Strikingly, at p, = 5.5po, 
maximum masses of the second and third family compact stars are nearly equal in the static limit. 
We find similar outcome for the other cases (not shown here). Thus, it seems there exist a critical 
value of Pt for a given A beyond which third family of compact stars with CFL core tend to 
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disappear with increase in pt. Below the critical value of p,, such hybrid stars in the non-rotating 
limit have maximum mass higher than their counterpart composed of hadronic matter. Our these 



56, 



results are substantiated by the earlier calculations performed using different EOSs [1 321 
The EOS used in Ref. [32] yields third family of compact stars in the static as well as in the 
mass-shedding limit sequences. For this EOS, the maximum mass of the compact star belonging 
to third family is larger by about 0.1 compared to its second family counterpart. On the other 



hand, for the EOSs used in Refs. [15 



57h . the third family compact stars exist in the static limit 



and disappears from the mass-shedding limit sequences. For these EOSs, maximum mass of the 



second and third family compact stars are nearly equal. In Ref. [|5qO, the maximum masses for the 



second and third family compact stars are found to be 1.57Mo and 1.55Mo, respectively. In Ref. 
[|57|], the maximum masses for the second and third family compact stars are 1.36Mo and 1.38Mo, 
respectively. 



C. Critical rotation frequency 

We have computed the values of the maximum or the critical rotation frequency /cnt for the 
stable configurations of hybrid stars with CFL core. The stable configurations of the compact stars 
rotating at a given frequency / satisfy. 

The Eq. ([231) is satisfied only for / < /cnt. To locate the critical frequency we first obtained the 
variation in the mass as a function of 6^ at fixed frequencies in steps of 50Hz. Then, for appropriate 
interval of the frequency, the calculations were repeated by varying the frequency in steps of 5Hz 
to determine the value of /cnt- In Figs. [51and[lUlwe plot the M - i?eq curves at fixed values of 
the rotational frequency. The black solid lines represent the results obtained at the / = /cnt. For 
the clarity, we mainly focus on the regions of the M - Rgg curves corresponding to the sequences 
of the hybrid stars which are relevant in the present context. In this region, the value of R^q 
decreases with increase in 6c. The results presented in Fig. l9l correspond to the cases for which 
third family of compact stars exists in the static as well as in the mass-shedding limit sequences. 
In Fig. [lUl we consider the cases for which third family of compact stars exist in the static limit, 
but, disappears from the mass-shedding limit sequences. We see that the value of /cnt, for the cases 
presented in Fig. l9l are larger than the highest observed rotation frequency 1 122 Hz. The value of 
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/crit - 775 - 900 Hz for the cases presented in Fig. \T0\ For p, = 4-po, the value of /cnt increases 
from 1370 - 1805 Hz as the CFL gap parameter A increases from 50 MeV to 150 MeV. It may be 



5711, but, the 



pointed out that the situation analogous to that of Fig. HOlis encountered in Refs. 11561 . 
value of /crit is about 350 - 650 Hz. 

In Table Hill we summarize the properties at the maximum mass of the hybrid stars with CFL 
core rotating with critical frequency. We compare the values of edfcnt) as given in this table with 
the maximum (6^(0)) and the minimum ( 6^(0)) energy densities at the centre of stable configura- 
tions of non-rotating hybrid stars. In the last two columns of Table Hill we give the values of dei 
and de2 calculated as, 

6,(0) 

and 

^'' = '- ^^^^ 

The values of 6,(0) and 6,(/crit) are taken from 6th and 4th columns of TablesHlandHlIl respectively. 
It is clear from the values of dei that the central energy density ec(/crit) is smaller than 6,(0) by 
25 - 30%. The values of de2 are noticeably larger for the cases considered in Fig. 13 than those of 
Fig. HOl For p, = 4po with A = 150 MeV, we get dei = 0.25 and de2 = 0.62. We would like to add 



that the EOS used in Ref. [I57l1 yields de2 - 0.4 for which third family of compact stars disappears 
from the mass-shedding limit sequence, but, exists in the static limit. It thus appears that de2 < 0.4 
disfavours the appearance of third family of compact stars in the mass- shedding limit sequence 
even if it exists in the static limit. It may be noted that the values of /cnt (Table Hill) is lower than 
the Kepler frequency (Table |lll) for the cases considered in Fig. |91 We would like to point out that 
the vlaue of /crit for a given EOS represents maximum rotation frequency for which the stability 
condition as given by Eq. (l23l) is satisfied. Whereas, the mass-shedding limit sequences can not 
be subjected to Eq. ([231) . Along these sequences, rotation frequency corresponds to the Kepler 
frequency which increases with the central energy density. This leads to the values of Kepler 
frequency at maximum mass higher than the fcrh for a given EOS. 



IV. CONCLUSIONS 



We construct static and mass-shedding limit sequences of hybrid stars for a set of EOSs ob- 
tained for diff"erent values of the CFL gap parameter and the deconfinement phase transition den- 
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sity. The hybrid stars considered are composed of CFL quark matter at the core, nuclear matter at 
the crust and mixed phase in the intermediate region. The hadronic part of the EOS is obtained 
using appropriately calibrated extended field theoretical based relativistic mean-field model. The 
EOSs of quark matter in the CFL phase corresponding to diff"erent values of the CFL gap param- 
eter and the deconfinement phase transition density are obtained using MIT bag model with an 
additional parameter that mimics the effect of including perturbative QCD corrections. The CFL 
gap parameter ranges from 50-150 MeV with the deconfinement phase transition density ranging 
from 4po - 6po (Po = 0.16 fm~^). 

We find the existence of the stable configurations of the static hybrid stars for all the differ- 
ent values of the CFL gap parameter considered with the deconfinement phase transition density 
4po - 5po. For the cases with CFL gap parameters 50 - 100 MeV with deconfinement phase tran- 
sition density 4po, the sequences of stable configurations of hybrid stars bifurcate from hadronic 
sequence when the central density exceeds the one at which the onset of mixed phase occur. In 
all the other cases, the stable configurations of hybrid stars form the third family of compact stars. 
When the deconfinement phase transition density is increased to 6po the stable configurations 
of hybrid stars exist only for CFL gap parameter A > 125 MeV. For the CFL gap parameter 
50 - 75 MeV with the deconfinement phase transition densities with 4po - 5po, the mass-radius 
relationship over broad range of mass and radius are in harmony with those deduced by applying 
improved hydrogen atmosphere model to fit the high quality spectra from compact star X7 in the 
globular cluster 47 Tucanae. The values of maximum mass 1.7 - 2.1 Mq for the hybrid stars are 
consistent with the currently measured maximum mass 1.76+0.20A/o of PSR J0437-4715 [|55l] . 

We find for several cases that the third family of compact stars disappear from the mass- 
shedding limit sequences, though, they appear in the corresponding static sequences. Our inves- 
tigation suggest that the third family compact stars is more likely to appear in the mass-shedding 
limit sequence provided they have maximum mass in the static limit higher than their counterpart 
composed of pure hadronic matter. Further, we have calculated the quantity de2 (Eq. [25l) obtained 
using the minimum and the maximum values of the central energy densities for the stable con- 
figurations of the static hybrid stars. The values of de2 is less than 0.4 for the cases for which 
third family of the compact stars disappears from the mass-shedding limit sequences, but, exists in 
the static limit. Except for these cases, the values of the critical rotation frequency for the hybrid 
stars with CFL core are larger than the highest observed frequency 1122 Hz. The relationship 
between the values of Sej and the disappearance of the third family of compact stars from the 
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mass-shedding limit sequences as observed in the present work is only preliminary. To establish 
this relationship, more investigations must be carried out using wide variety of EOSs. 
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TABLE I: The maximum mass of hybrid stars with CFL core and hadron stars in the static limit and cor- 
responding central energy density and radius obtained for different values of the CFL gap parameter A and 
the deconfinement phase transition density pt. 









Hadron Stars 




Hybrid Stars 




Pt 


A 


€ 


M 


R 


€ 


M R 




(Mev) (lO'^g/cm^) 


Mq 


(km) 


(lO^^g/cm^) 


Mo (km) 




50 


0.721 


1.63 


13.49 


2.716 


2.12 10.57 




100 


0.642 


1.45 


13.52 


2.863 


2.05 10.16 




150 


0.512 


1.04 


13.48 


3.244 


1.97 9.42 




50 


0.873 


1.94 


13.34 


2.863 


1.96 10.53 


5po 


100 


0.808 


1.83 


13.42 


3.072 


1.89 10.10 




150 


0.690 


1.62 


13.55 


3.663 


1.80 9.28 




50 


1.032 


2.10 


13.15 






6po 


100 


0.959 


2.04 


13.24 








150 


0.850 


1.88 


13.38 


4.071 


1.69 9.30 
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TABLE II: The central energy density, radius and Kepler frequency at the maximum mass for the hybrid 
stars with CFL core and hadron stars obtained for different values of the CFL gap parameter A and the 
deconfinement phase transition density p,. 











Hadron Stars 






Hybrid Stars 




Pt 


A 


e 


M 


R 




e 


M 


R iK 




(Mev) (lO'^g/cm^) 


Mo 


(km) 


(Hz) (10i5g/cm^) 


Mo 


(km) (Hz) 




50 


0.721 


2.09 


18.84 


1028 


2.383 


2.48 


13.95 1698 


4po 


100 


0.642 


1.84 


19.06 


953 


2.665 


2.38 


13.21 1803 




150 


0.512 


1.32 


19.24 


801 


2.934 


2.29 


12.18 1993 




50 


0.873 


2.45 


18.25 


1158 








5po 


100 


0.808 


2.32 


18.51 


1 107 


2.842 


2.17 


13.44 1687 




150 


0.690 


2.08 


18.95 


1019 


3.339 


2.05 


12.20 1889 



TABLE III: Properties at the maximum mass of hybrid stars with CFL core rotating at the critical frequency 
as obtained for different values of the CFL gap parameter A and the deconfinement phase transition density 
Pt- The values of de^ and de2 are obtained using Eqs. (I24l) and (l25l) . 



Pt 


A 


/crit 


e 


M 


/?eq dei de2 




(Mev) (Hz) (lO^^g/cm^) (Mo) (km) 




50 


775 


2.160 


2.01 


11.58 0.25 0.35 


5po 


100 


1125 


2.205 


1.98 


11.49 0.28 0.45 




150 


1540 


2.685 


1.94 


11.30 0.27 0.55 


6po 


150 


905 


2.951 


1.74 


10.62 0.28 0.30 
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FIG. 1 : (Color online) Pressure as a function of energy density for different values of the CFL gap parameter 
with hadron to CFL quark matter phase transition densities = 4po (left), 5po (middle) and 6po (right)- 
The EOS for pure hadronic matter is shown by the black dotted hne. 
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FIG. 2: (Color online) Variations of bag constant as a function of the CFL gap parameter at fixed values of 
hadron to CFL quark matter phase transition densities 4po - 6po- 
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FIG. 3: (Color online) Variations of the average value of quark chemical potential at the phase transition 
density as a function of the CFL gap parameter. 
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FIG. 4: (Color online) Plots for the mass-radius relationships for the static sequences obtained using fixed 
values of the CFL gap parameter A ranging from 50-150 MeV with deconfinement phase transition density 
pt - 4po- The dotted black line represent the static sequence of compact stars composed of the hadronic 
matter. The solid circle on each of the curves denote the end of the mixed phase. The curves on the left 
to the solid circles represent the hybrid stars with CFL core. The dot-dashed maroon curves are the mass- 
radius contours deduced with 90% confidence by fitting the high quahty spectra from the compact star X7 



in the globular cluster 47 Tucanae i49h . 
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FIG. 5: (Color online) Same as Fig. |4l but, for the deconfinement phase transition density p, = 5po. 
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R [km] 

FIG. 6: (Color online) Same as Fig. HI but, for the deconfinement phase transition density p, = 6po. Inset 
highlights the appearance of third families of compact stars for the CFL gap parameter A > 125 MeV. 
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FIG. 7: (Color online) Relationship between mass M and the circumferential equatorial radius /?eq for mass- 
shedding limit sequences for different values of the CFL gap parameter and the phase transition densities 
4po (upper panel), 5po (middle panel) and 6po (lower panel). The curves on the left to the solid circles 
represent the hybrid stars with CFL core. The symbol cross (x) on the different curves marks the maximum 
mass of the hybrid star. 
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FIG. 8: (Color online) Plots for the mass-shedding limit sequences (upper panel) and the static sequences 
(lower panel) for the CFL gap parameter A = 150 MeV with deconfinement phase transition density = 
5po-6po. The parts of the curves left to the solid circles in the upper and lower panel represent the sequences 
of hybrid stars with the CFL quark matter core. 
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FIG. 9: (Color online) Plots for the mass verses circumferential equatorial radius /?eq at fixed values of the 
rotational frequency. The black solid lines represent the results obtained at the critical frequencies fcnt- For 
/ > /crit> third families of compact stars do not exist. The values of pt and A considered are such that they 
yield third family compact stars in the static as well as in the mass-shedding limit sequences. 
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